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ABSTRACT 

We have studied the rapid X-ray time variabiUty in 149 pointed observations with the Rossi X-ray 
Timing Explorer (RXTE)'s Proportional Counter Array of the atoU source 4U 1636-53 in the banana 
state and, for the first time with RXTE, in the island state. We compare the frequencies of the 
variability components of 4U 1636-53 with those in other atoll and Z-sources and find that 4U 1636- 
53 follows the universal scheme of correlations previously found for other atoll sources at (sometimes 
much) lower luminosities. Our results on the hectohertz QPO suggest that the mechanism that sets 
its frequency differs from that for the other components, while the amplitude setting mechanism 
is common. A previously proposed interpretation of the narrow low-frequency QPO frequencies in 
different sources in terms of harmonic mode switching is not supported by our data, nor by some 
previous data on other sources and the frequency range that this QPO covers is found not to be 
related to spin, angular momentum or luminosity. 

Subject headings: accretion, accretion disks — binaries: close — stars: individual (4U 1636-53,411 
1820~30,4U 1608-52,4U 06144-09,4U 1728-34) — stars: neutron — X-rays: stars 
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1. INTRODUCTION 

Low-mass X-ray binaries (LMXBs) can be divided 
into systems containing a black hole candidate (BHC) 
and those containing a neutron star (NS). The accre- 
tion process onto these compact objects can be stud- 
ied through the timing properties of the associated X- 
ray emission (see, e.g., ly an der Klis 2006, for a review). 
iHasinger fc van der Klis[(| 19891 ) classified the NS LMXBs 
based on the correlated variations of the X-ray spec- 
tral and rapid X-ray variability properties. They dis- 
tinguished two sub-types of NS LMXBs, the Z sources 
and the atoll sources, whose names were inspired by the 
shapes of the tracks that they trace out in an X-ray color- 
color diagram on time scales of hours to days. The Z 
sources are the most luminous; the ato ll sources cover a 
much wider range in luminosities (e.g. . iFord et al.ll2000l . 
and references therein). For each type of source, several 
spectral/timing states are identified which are thought 
to arise from qualitatively different inner flow configura- 
tions. In the case of atoll sources, the main three states 
are the extreme island state (EIS), the island state (IS) 
and the banana branch, the latter subdivided into lower- 
left banana (LLB), lower banana (LB) and upper banana 
(UB) states. Each state is characterized by a unique 
combination of color color diagram and timing behav- 
ior. The EIS and the IS occupy the spectrally harder 
parts of the color color diagram (CD) corresponding to 
lower X-ray luminosity {Lx). The different patterns they 
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show in the CD are traced out in days to weeks. The 
hardest and lowest state is generally the EIS, which 
shows strong low-frequency flat-topped noise. The IS 
is spectrally softer than the EIS. Its power spectrum is 
characterized by broad features and a dominant band- 
limited noise (BLN) component which becomes stronger 
and lower in characteristic frequency as the flux decreases 
and the > 6 keV spectrum gets harder. In order of in- 
creasing Lx we encounter the LLB, where the twin kHz 
QPOs are first observed, the LB, where dominant 10-Hz 
BLN occurs and finally, the UB, where the < 1 Hz (power 
law) very low frequency noise (VLFN) dominates. In the 
banana states, some of the broad features observed in 
the EIS and the IS become narrower (peaked) and oc- 
cur at higher frequency. The twin kHz QPOs can be 
found in LLB at frequencies in excess of 1000 Hz, only 
one is seen in the LB, a nd no kHz QPOs are detected 
in the UB ( see reviews by Hasin ger fc van der Klislfl989l : 
Ivan der Kl is 2000, 2004, 2006, for detailed descriptions 
of the different states) . 

4U 1636-53 is an atoll source (jHasinger fc van der KlisI 
I1989D which has an or bital period of ^ 3.8 hours 
(jvan Paradijs et ahl Il990( ) and a companion star with 
a i nass of ^ 0.4 M c^ (assuming a NS of ~ 1.4 Mq, 
see iGiles et al] l2002l for a discussion) . It was first ob- 
served as a strong co ntinuous X-ray sourc e (Norma X- 
1) with Copernicus (jWillmore et all Il974[ ) and Uhuru 
(jCiacc oni et al. lll974D. 4U 1636-53 is an X-ray burst 
source ( Hoff man et al.l I1977D which shows asymptotic 
burst oscillation frequencies of ~ 581 Hz (see e.g. 
IZhang et al]ll997t iGiles et al.ll2002[ ). Th is is probably 
the approximate spin frequency; although iMilleil ()1999D 
presented evidence that these oscillations might actually 
be the second harmonic of a neutron star spin frequency 
of ^ 290 Hz , this w a s not confirmed in fur t her w ork by 
IStrohmaved (|2001a[ ). iPrins fc van der KlisI (Il997l) stud- 
ied the aperiodic timing behavior of 4U 1636-53 with 
the EXOSAT Medium Energy instrument up to frequen- 
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cies of ^ 100 Hz both in th e island and the banana 
state. IWiinands et al.l (|1997| ). using observations with 
RXTE, discovered two simultaneous quasi-periodic os- 
cillations (QPOs) near 900 Hz and 1176 Hz when the 
source was in the banana state. The frequency difference 
IS.V between the two kHz QPO peaks is nearly equal to 
half the burst oscillation frequency, similar to what has 
been observed in other sources with burst oscillations 
or pulsation frequency > 400 Hz. To the extent that 
this implies Az/ ~ Vspin/'^, thi s is inconsistent with spin- 
orbit beat-frequency models (IWiinands et al.i,2003') fo r 
the kHz QPOs such as proposed bv iMiller et al 
Other complications for beat frequency models include 
the fact that At/ is neither constant (e.g. in Sco X- 
1, Ivan der Klis et all Il997f ) nor exactly equal to half 
the burst oscillation frequency. Generally, Ai/ decreases 
as the kHz QPO frequency increases, and in 4U 1636- 
53, observations have shown Ai/ at frequencies lower as 
well as higher than half the burst oscilla tion frequency 
(jMendez et al." 1998a; 'Jonker et al."2002a'). 

[van Straaten ct al. (2002, 2003 ) compared the timing 
properties of 4U 0614-^09, 4U 1608-52 and 4U 1728- 
34 and conclude that the frequencies of the variability 
components in these sources follow the same pattern of 
correlations. ' Di Salvo et al.l (|2003f ). based on five detec- 
tions of kHz QPOs in 4U 1636-53 in the banana state was 
able to show that at least in that state the source might 
fit in with that same scheme of correlations. The de- 
tailed investigation of 4U 1636-53 is import a,nt because 
it is one of the most luminous atoll sources ()Ford et al.l 
l2000t ) that shows the full complement of island (this pa- 
per) and banana states and that also shares other timing 
feature s with often less l umino us atoll sources. For ex- 
ample, iRevnivtsev et al.l ()2001[ ) found a new class of low 
frequency QPOs in the mHz range which they suggested 
to be associat ed with n uclear burning in 4U 1636 53 and 
4U 1608-52. iMended (j2000l ) and Mcnd ez et all ((20011) 
compared the relations between kHz QPOs and inferred 
mass accretion rate in 4U 1728-34, 4U 1608-52, Aql X- 
1 and 4U 1636-53, and showed that the dependence of 
the frequency of one of the kHz QPOs upon X-ray inten- 
sity is c omplex, but similar among sources. iJonker et al.l 
(jlOOOi) discovered a third kHz QPO in 4U 1608-52, 
4U 1728-34, and 4U 1636-53 which is lik ely an upper 
side band to the lower kHz QPO. Recently, iJonker et al.l 
(1M)5) found in 4U 1636-53 an additional (fourth) kHz 
QPO, likely the corresponding lower sideband. 

In this paper, we present new results for low frequency 
noise with characteristic frequencies 1—100 Hz and QPOs 
in the range 100 — 1260 Hz, for the first time including 
RXTE observations of the island state of this source. 
These results better constrain the timing behavior in 
the various states of 4U 1636-53. We compare our re- 
sults mainly with those of the atoll sources 4U 0614-1-09, 
4U 1608-52 and 4U 1728-34 and find that the frequency 
of the hectohertz component may not be constant as 
previously stated, but may have a sinusoidal like mod- 
ulation within its range from ~ 100 to ~ 250 Hz. Our 
results also suggest that the mechanism that sets the fre- 
quency of the hHz QPOs differs from that for the other 
components, while the amplitude setting mechanism is 
common. Finally, we demonstrate that it is not possi- 
ble to clearly distinguish between two harmonics of the 
low-frequency QPO L^p across different sources, as was 



previously thought ijvan Straaten et"al]|2003f) . 

2. OBSERVATIONS AND DATA ANALYSIS 

We use data from the Rossi X-ray Timing Explorer 
(RXTE) Proportional Counter Array (P CA; for instru- 
ment information see IZhang et al.lll993f ). There were 
149 pointed observations in the four data sets we used 
(60032-01, 60032-05, 70036-01, 80425-01 & 90409-01), 
each consisting of a fraction of one to several entire satel- 
lite orbits, for ~ 1 to ^ 26 ksec of useful data per 
observation. We use the 16-s time-resolution Standard 
2 mode data to c alcula te X-ray colors as described in 
lAltamirano et all (|2005( ) . Hard and soft color are defined 
as the 9.7-16.0 keV / 6.0-9.7 keV and 3.5-6.0 keV / 2.0- 
3.5 keV count rate ratio, respectively, and intensity as the 
2.0-16.0 keV count rate. Type I X-ray bursts were re- 
moved, background was subtracted and deadtime correc- 
tions were made. In order to correct for the gain changes 
as well as the differences in effective area between the 
PCUs themselves, we normalize o ur colors by the cor- 
responding Crab color v alues (see iKuulkers et al.l 11994 
Ivan Straaten et 1111120031 ) that are closest in time but in 
the same RXTE gain epo ch, i.e., with the sam e high volt- 
age setting of the PCUs (jJahoda et al.ll2005l ). In Table [2] 
we show for reference the all-epoch averaged colors for 
Crab Nebula. All active PCUs were used to calculate 
the colors in 4U 1636-53 except for observation 60032- 
01-01-02, where due to a PCU3 malfunction we only used 
PCUs and 2. Figure [T] shows the color-color diagram of 
the 149 different observations that we used for this anal- 
ysis, and Figure |2| the corresponding hardness-intensity 
diagrams (soft and hard color vs. intensity). 

For the Fourier timing analysis we used data from 
the ~ 125/xs (1/8192 s) time resolution Event mode 
E_125us_64M_0_ls. First, we used a 2 second-binned 
light curve in order to detect and remove data drop-outs 
and X-ray bursts (these data were also excluded from 
the rest of the analysis). Leahy-normalized power spec- 
tra were constructed using data segments of 128 seconds 
and l/8192s time bins such that the lowest available fre- 
quency is 1/128 w 8x lO^'' Hz and the Nyquist frequency 
4096 Hz. No background or deadtime corrections were 
made prior to the calculation of the power spectra. We 
first averaged the power spectra per observation. We in- 
spected the shape of the average power spectra at high 
frequency ( > 2000 Hz) for unusual features in addition 
to the usual Poisson noise. None were found. We then 
subtracted a Poisson noise spectrum estimated from the 
power between 3000 and 4000 Hz, where neither intrin- 
sic noise nor QPOs are known to be present, using the 
method developed by Klein- Wolt (2004) based on the 
analytical function of Zhang et al.l (|l995f l. The resulting 

S lower spectra were converted to squared fractional rms 
van der KlisI |1995[ ). In this normalization the square 
root of the integrated power density equals the variance 
of the intrinsic variability in the source count rate. In 
order to improve the statistics, observations were aver- 
aged together if they described the same source state. 
Since it is known from previous work on similar sources 
that the position of the source in the color-color diagram 
generall y is well correlated to its spectral/timing state 
(see e.g. Ivan de r Klis 2006, and references within), we 
first grouped observations with similar colors. Within 
each group, we then compared the shape of each aver- 
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Fig. 1. — Hard color versus soft color normalized to the Crab Nebula as explained in Section [T] Each circle represents the average 
soft/hard color of one of the observations used for this paper. The filled triangles mark averages of 1 to 32 observations and are labeled 
with letters, in order from {A) the island state, Lower Left Banana (J) to the Lower banana (N). 
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Fig. 2. — Soft color vs. intensity (left) and hard color vs. intensity (right) in units of the Crab Nebula as explained in Section[T] Symbols 
as in Figure [T] For clarity, the dashed line separates the observations corresponding to the IS (left), from the observations corresponding 
to the BS (right). 



age power spectrum with all of the other ones to cre- 
ate subgroups in which all power spectra had a depen- 
dence of power on frequency that was identical within 
errors. So, narrow features had to be at the same fre- 
quency for average power spectra to be added together. 
The resulting data selections are labeled interval A to 
N (see Table [T] for details on which observations were 
used for each interval and their colors). A disadvan- 
tage of this method is that we can loose information 
about narrow features moving on time scales shorter than 
an observation, such as the lower kilohertz QPO (see 



e.g. iBerger e~a l. 1996; Di Salvo et al.ll2003D . The "shift 
and add" method (Mendez et al. l ll998bD~ to some extent 
might be able to compensate for this; we explore in the 
Appendix this issue. Our method is the best suited one 
to study the behavior of the broad features such as typi- 
cally seen in low mass X- r ay binaries' power spectra (e.g. 



van Straa ten et all f200l l2003l . 120051 : lAltamirano et all 
2005; Lin ares et al.ll2005h . For these broad components, 
which are the main aim of this paper, the gain in signal 
to noise due to this averaging process outweighs a minor 
additional broadening due to frequency variations. 
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Interval Al 


Observation 


Soft color (Crab) 


Hard color (Crab) 


Intensity (Crab) 


80425-01-04-01 


1.2306 ±0.0034 


1.0493 ±0.0032 


0.0557 ± 0.0001 


Interval A2 


80425-01-03-00 
90409-01-01-00 
90409-01-01-01 
90409-01-02-00 


1.2097 ±0.0056 
1.1945 ±0.0049 
1.1954 ±0.0032 
1.2002 ±0.0040 


1.0585 ± 0.0053 
1.0590 ± 0.0049 
1.0403 ± 0.0031 
1.0424 ± 0.0042 


0.0400 ± 0.0001 
0.0370 ± 0.0001 
0.0418 ± 0.0001 
0.0487 ± 0.0001 



TABLE 1 

Observations used for the timing analysis. The colors and intensity are normalized to Crab (See Section [2J. The complete 

TABLE CAN BE OBTAINED DIGITALLY FROM ApJ. 



PCU Number 


Crab's soft color 


Crab's hard color 


Crab's intensity (c/s) 





2.21 ± 0.04 


0.56 ±0.009 


2552 ± 21 


1 


2.27 ±0.03 


0.55 ±0.008 


2438 ±21 


2 


2.22 ± 0.04 


0.58 ±0.010 


2424 ± 20 


3 


2.42 ± 0.03 


0.57 ±0.009 


2365 ± 20 


4 


2.34 ±0.03 


0.58 ±0.011 


2299 ± 21 



TABLE 2 

Average soft, hard and intensity of the Crab Nebula over all epochs and per PCU. Note that we have used averaged 

VALUES per day IN OUR ANALYSIS; THE NUMBERS LISTED HERE ARE REPRESENTATIVE FOR THOSE DAILY AVERAGES. We QUOTED THE 
averaged QUADRATIC ERRORS, I.E. \/Y17 ^^i/™, WHERE Ak IS EITHER THE SOFT, HARD OR THE INTENSITY DAY ERROR. 



To fit the power spectra, we used a multi-Lorentzian 
function: the sum of several Lorentzian components plus, 
if necessary, a power law to fit the very low frequency 
noise at < 1 Hz. Each Lorentzian component is denoted 
as Li, where i determines the type of component. The 
characteristic frequency {vmax as defined below) of Li is 
denoted I'i. For example, L„ identifies the upper kHz 
QPO and its characteristic frequency. By analogy, 
other components have names such as Lg (lower kHz), 
LhHz (hectohertz), Lh (hump). Lb (break frequency), 
and their frequencies are ve, v^hz, Vh and Vh, respec- 
tively. For reference, in Figure [3] we show two repre- 
sentative power spectra in which we labeled the differ- 
ent components. Using this multi-Lorentzian function 
makes it straightforward to directly compare the differ- 
ent components in 4U 1636-53 to tho se in previous works 
which used the same fit function fe.g..|Belloni et al1l2002l: 
van Straaten et all l200i I2003L l2005i : lAltamirano et all 
20051 and references therein). 

We only include those Lorentzians in the fits whose 
single trial significance exceeds Scr based on the nega- 
tive error bar in the power integrated from to cx) (i.e. 
we include only those Lorentzians whose integral power 
is at least 3 times higher than zero based on the neg- 
ative Icr error) and whose inclusion gives a > 3(t im- 
provement of the fit according to an F-test. We give 
the frequency of the Lorentzians in ter ms of charac- 
teristic frequency Vmax as introduced by iBelloni et al.l 
(EOOl): 

For the quality factor Q we use the standard definition 
Q = vq/FWHM. FWHM is the full width at half max- 
imum and vq the centroid frequency of the Lorentzian. 
Note that Q values in excess of ~ 3 will generally be af- 
fected by smearing in an analysis such as ours. Such val- 
ues are commonly seen in Llf, Li and Lu- In Section [3] 
we indicate in which cases this could have occurred. 

We only report the results for Vmax ^ 1 Hz. 4U 1636- 
53 is one of three atoll sources which are known to show 
milihertz QPOs whi ch affect the power law behavior of 
the noise at < 1 Hz (|Revnivtsev et aLll2001|) . A different 
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Fig. 3. — Representative power spectra of the island state (above 
- interval B) and the banana state (below - interval J) with their 
components. 

kind of analysis is needed to study th ese QPOs; we will 
repor t the results in a separate paper (jAltamirano et aU 
[2008I) . 

3. RESULTS 

Figures [1] and [5] show that in order A to H, the spec- 
trum becomes softer (i.e. hard and soft color both de- 
crease), and the intensity changes little. From H to L 
the soft color remains approximately constant but above 
6 keV the spectrum becomes even steeper (i.e., the hard 
color decreases further) and, from interval G, the inten- 
sity increases. Finally from L to N, below 6 keV the spec- 
trum becomes flatter and above 6 keV it remains approx- 
imately constant in slope, while the intensity continues 
increasing. Similar behavior has been observed in other 
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Interval 




FWHMt,,^ 






Auo 


Ratio Vuo/vtd 


H 


860.4 ± 1.7 


136.2 ±4.8 


565.4 ±5.1 


90.4 ± 12.3 


295.0 ±5.4 


1.52 ±0.01 


I 


896.8 ±2.6 


115.8 ±6.3 


580.4 ± 5.1 


77.2 ± 12.5 


316.4 ±5.7 


1.54 ±0.01 


J 


972.1 ±3.5 


98.0 ±8.7 


661.1 ± 1.8 


73.8 ±4.1 


311.0 ±3.9 


1.470 ±0.006 


K 


992.9 ±6.6 


177.4 ± 16.9 


718.3 ± 2.4 


121.5 ± 6.3 


274.6 ±5.7 


1.38 ±0.01 


L 


1147.2 ± 16.6 


123.4 ±27.8 


836.1 ±0.7 


64.1 ± 1.9 


311.1 ± 16.6 


1.37 ±0.02 



TABLE 3 

Central frequencies, full width at half maximum (FWHM) and the frequency difference Ai^q between the kHz QPOs for 

THE 5 intervals WHERE BOTH kHz QPOS WHERE DETECTED SIGNIFICANTLY (> 3(t). *: NOTE THAT THESE VALUES MIGHT HAVE BEEN 

affected by SMEARING (SEE TEXT) . 



atoll sources which are moving from the island to the 
lower left ba nana and then to the lowe r banana state (see 
for e xample Ivan Straaten et aTl 120031 : lAltamirano et ahl 
[2005h . 

In Figure 2] we show the average power spectra A — N 
with their fits. Two to five Lorentzian components were 
needed for a good fit, except in power spectrum I, where 
an extra component is needed, for a total of six. Ta- 
ble [5] gives the fit results. Power spectra Al and A2 
have the same within errors, and only slightly differ- 
ent colors. We treat these two power spectra separately 
because Al could be fitted with 5 significant components 
and A2, as well as the combined spectrum A1+A2, only 
with 3. Note, that power spectrum Al is the average of 
one observation (See Table[T|) which was performed in be- 
tween the observations used for power spectrum A2. In 
Figure we show our measured characteristic frequency 
correlations (black) together with those previously mea- 
sured in other atoll sources (grey). In intervals H-L, the 
twin kilohertz QPOs (L„ and L^) are identified unam- 
biguously. The correlation between the lower and the 
upper kilohertz QPO is the same as that found in th e 
other atoll sources studied bv lvan Straaten et al.l (|2003D . 
For intervals M and N, only L„ is observed. 

We separately measured the centroid frequencies vq of 
the kHz QPOs (see Table [3]). The centroid frequency 
difference /S.v varied between 274.6 ± 5.7 Hz (in interval 
K) and 316.4 ± 5.7 Hz (inte rval I). These values are be- 
tween the extre mes found by l Jonker et al.l ()2002aD (Az/ = 
330±9 Hz) and lDi Salvo"et a l. (2003) (Az/ = 242±4 Hz). 
Note that those authors used much shorter time intervals 
to detect the kHz QPOs and hence were more sensitive 
to short lived extreme cases. Our average power spectra 
contain those data which is necessary to detect the broad 
components well, and hence our measured Ai/'s occur at 
intermediate value s, which is when th e kHz QPOs are 
strongest (see e.g. iMendez et ahllMl . So our method 
averages out the extreme cases. In power spectrum N, L„ 
reaches the highest centroid frequency found among the 
intervals analyzed, 1259 ± 10 Hz. Note that power spec- 
tra M and N are the result of averaging large amounts 
of data with no significant kHz QPOs in individual ob- 
servations 2.2 X lO"' seconds and ~ 2.5 x 10^ seconds, 
respectively), based on the position of 4U 1636-53 in the 
color-color diagram. Figure [5] displays the upper kHz 
QPOs in power spectra M and N more clearly. 

Between ~ 100 Hz and ^ 200 Hz, a Lorentzian with 
quality factor Q ~ 1 is often found in atoll sources (see 
Ivan der Klisll200^ . and references therein). This feature 
is called hectohertz QPO and, in contrast to the other 
components, its frequency remains confined to this rela- 
tively narrow range as increases. When Vu ^ 800 Hz, 
the twin kilohertz QPOs are usually identified unambigu- 



ously, and so is the hectohertz QPO. For ^ 600 Hz, 
the lower kHz QPO could also have frequencies between 
~ 100 Hz and ^ 300 Hz if it would be present, which 
makes it difficult to classify the QPOs fovmd in that range 
as either Li, L^hz or a blend without more information. 
In our data, this is the case for intervals A to E; in Table[5] 
and hereafter we identify those Lorentzians as hectohertz 
QPOs. This identification is supported by the fact that 
in intervals F and G, i.e., for 600 < f„ < 800 Hz, 
is undetected; this component seems to appear only at 
i^u > 800 Hz. Interval I shows a ~ S.ltr (single trial) 
peak with lymax = 229 ±9 Hz which is a factor ~2 higher 
than the usual hHz in that range. This QPO may be the 
second harmonic of the hectohertz QPO simultaneously 
found at a characteristic frequency of 113 ±4 Hz (see Fig- 
ure[7l however, note that this feature must be interpreted 
with care due to its low -3.1tT- single trial significance). 
As a result of refitting the power spectrum using centroid 
frequencies, we find that the second harmonic QPO is at 
vq = 228 ± 10 Hz while the first harmonic hHz QPO is 
at 1^0 = 107 ± 5 Hz for a ratio of 2.13 ± 0.13, consistent 
with 2. 

Figure [5] shows that also Lh and Lf, lie on the corre- 
lations previously observed in other sources. However, 
our results show that i>b may anti-correlate with at 
Vu ^ 1100 Hz. This result has already been observed in 
Z sources, however for atoll sources this behavior has not 
been observed with certainty (See Section 2]). Lb2 seems 
to have lower frequencies than in other atoll sources. To 
further investigate this, in Figure [8] we plot vi,2 versus 
Vu with different symbols for each of the 4 atoll sources 
for which this component has been measured. Clearly, 
the range in which has been found for similar i/u 
is rather large (up to nearly a decade), particularly at 
Vu < 1000 Hz. We also studied the possibility that the 
rms of Li,2 could be related with its frequency, but no 
relation was found. 

Intervals A, B, D and E show a narrow QPO 
with a characteristic frequency between Vb and Vh 
(see Table For other neutron st ars such nar- 

row QPOs were previously rep orted by 'Yoshida et al.l 
(1993i) in 4U 1608-52, by Bc Uoni et ai.l (2002) in the 
low lumi nosity bursters IE 1724-3 045 and OS 1826- 
24, by Ivan Straaten et al.l (l2002t) in 4 U 0614+09, 
4U 172 8-34. bvlvan Straaten et all pOOSi ) in 4U 1608- 



52, by lAltamirano et a" 



Ivan Straaten et all (|2005l ) 



( 2005ir m 
in the 



4U 1820-30, by 
accreting millisec- 
ond pulsars (AMP) XTE J0929- 314, XTE J18 1 4-338 
and SAX J1808.4-3658 and by ILinares et all (pOOSh 
in XTE J1807-294. Simila r features were al s o seen 
in the BHCs Cyg X-1 by iPo ttschmidLefap ()2003f ) 
and GX 339-4 by iBelloni etal. ( 2002„ but also see 
Ivan Straaten et al. 1 120031) . Following Ivan Straaten et ahl 
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Fig. 4. — Power spectra and fit functions in the power spectral density times frequency representation. Each plot corresponds to a different 
region in the color-color and color-intensity diagrams (See Figures [T] and |2]l . The curves mark the individual Lorentzian components of the 
fit. For a detailed identification, see Table [5l and Figure |5] 




Vu (Hz) 

Fig. 5. — The characteristic frequencies Umax of the various power spectral components plott ed versus Uu- The black bullets mark 
the results for 4U 1636-53. The other symbols mark the atoll sources 4U 0614+09, 4U 1728-34 llvan Straaten et al.ll200g) . 4U 1608-52 
jraTn Straatcn ct al. 2003), and Aql X-1 (Rcig ct al. 2004) and the low luminosity bursters IE 1724-3045, GS 1826-24 and SLX 1735-269 
Jvan Straaten ct al...2005i . but also see .BcUoni ct al...2002V 
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Fig. 6. — kHz QPOs of intervals M and A'', respectively. These 
are Leahy normalized power spectra with no Poisson noise subtrac- 
tion. 



'lli ' ' ^ 

T hHz QPO Twin kHz QPOs 




Fig. 7. — Part of the power spectrum of interval I, showing the 
twin kHz QPOs, the hHz QPO and its possible harmonic. This 
is a Leahy normalized power spectrum with no Poisson noise sub- 
traction. 

(|2003l ). for clarity we have omitted these QPOs (^lf) 
from Figure [5] In Figure [H] we plot their character- 
istic frequencies vs. i^h- The results for 4U 1636-53 
are in the range of, but seem to follow a different re- 
lation than, the two rela tions previously suggested by 
Ivan Straaten et al] (|2003l) based on other sources. The 
Llf QPOs in 4U 1636-53 cannot be significantly de- 
tected on timescales shorter than the duration of an av- 
erage observation. 

For completeness, in Figure \W\ we plot both the frac- 
tional rms amplitude and quality factor of the L^p com- 
ponent versus i>u ■ Although the fractional rms amplitude 
of 4U 1636-53 increases with v^, no general trend is ob- 




4U 0614+09 
4U 1728-34 
4U 1608-52 V 

i [ ^ ^ ^ ^ 4U 1636-53 ^^T^ " 

500 600 700 800 900 1000 1100 1200 1300 
v„ (Hz) 

Fig. 8. — ui,2 v ersus Uu for the 4 atoll sources 
4U 0614-1-09 , 4U 1728-34 llvan Straaten et al .lf200g). 4U 1608-52 
||van Straate n et al.|[2003l) and 4U 1636—53 (this paper). Note that 
the ope n square at Vu ~ 1233 Hz could be interpreted as either L), 
or Li,2 l lvan Straaten et al.ll200^ . 

served among the 7 sources shown in this Figure. The 
quality factor Qlf (which may have been affected by 
smearing, see Section [2]) seems to be unrelated to Vu for 
all the sources shown. 

In Figure [TT] we plot the fractional rms amplitude of 
all components (except L^p, see Figure [TO]) versus 
for the atoU sources 4U 0614-^09, 4U 1728-34, 4U 1608- 

52 and 4U 1636-53. The rms of the upper kHz QPO 
for all sources approximately follows the same trend: it 
increases up to I'u ~ 750 — 800 Hz, and then starts to de- 
crease. This seems also to be the case for LhHz and Li,. 
Except for 4U 0614-f 09, the data suggest that at i^u ?J 
1100 Hz the rms of Lu does not decrease further, but 
remains approximately constant. At ~ 750 — 800 Hz, 
the rms amplitudes of L^hz and Lb start to decrease 
(see also Ivan Straaten et aI]l2003D . The rms of Lh of 
4U 1636-53 also seems to follow the general trends ob- 
served for the other atoll sou rces. Some of these result s 
were previo usly reported by van Straaten et al.' (2003'), 
iMendez et al. (2001) and Barret et al. (2005a). 4U 1636- 

53 stands out by the fact that the rms of Lf,2 and LhHz at 
Vu < 900 Hz is always smaller than in the other sources. 
Moreover, and again contrary to the other sources, in 
4U 1636-53 the rms of Lh2 remains approximately con- 
stant as Vu increases from ~ 800 Hz to ~ 1200 Hz. 

In Figure [T2| we plot the quality factors Q„, and 
QhHz versus f„. As noted in Section |2l the Q values 
of Lg, and to a lesser extent, L„, have likely been af- 
fected by smearing. The Q values of the other compo- 
nents are not plotted since they are usually broad. The 
data on 4U 1636-53 are in general agreement with what 
was found using similar methods on the other sources 
(|van Straaten et aLll2002l ). Qu increases monotonically 
with Vu until Vu ~ 900 — 1000 Hz. At this frequency, 
Qu seems to decrease for all sources, to immediately in- 
crease again as Vu increases. Q j shows a rather ra ndom 
behavior due to smearing (see lDi Salvo et al.ll2003 and 
iBarret et al.l l2005al |bl for Qi measurements less affected 
by smearing) ; we display this quantity in Figure [T^ for 
comparison to previous works using the same method. 

QhHz shows a complicated behavior. To further inves- 
tigate this behavior, we first re-binned the Qhhz data by 
a factor 3 and fitted a straight line. The best fit gives 
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10 100 
Vh (Hz) 



Fig. 9. — (Above) Characteristic frequencies ui^p and V]^pi2 
(see text) versus Vf^ . The symbols are labeled in the plot, and rep- 
resent the frequencies of the QPOs in the atoll sources 4U 1728- 
34, 4U 0614+09, 4U 1608-52, 4U 1820-30, the BHCs Cyg X-1 
and GX 339-4, the low luminosity bursters IE 1724-3045 and 
GS 1826-24 and the accreting millisecon d pulsars XTE J0929-314 
XTE J1814-338 and SAX J1808. 4-3658 l lvan Straaten et al.ll2003l . 
[200 5: Altamirano ot al. 2005). The filled triangles show the results 
for 4U 1636-53. The drawn line indicates a power-law fit to the 
vlf vs. Vh relation of the low-luminosity bursters IE 1724—3045 
and GS 1826-24, and the BHC GX 339-4. The dashed line is 
a power law with a normalization half of that of the drawn line. 
(Bottom) A zoom of the high frequency region. 

/dof = 93/12 ~ 7.7. Since the data appears to show 
two bumps separated by a minimum creating a roughly 
sinusoidal pattern, we also tried to fit a straight line plus 
a sine wave. The best fit has a x^/dof = 16.9/9 ~ 1.8 
which gives a 3.4(t improvement of the fit based on 
an F-test. In Figure O it can be seen that VhHz ap- 
pears to display a similar pattern. Fitting the rela- 
tion of VhHz versus Vu with only a straight line gives 
a /dof = 156/12 — 13 while a straight line plus a sine 
gives a x^/dof = 28.4/9 ^ 3.15, once again, we find a 
3.4(7 improvement of the fit based on an F-test. The re- 
sults of the sinusoidal fits, in which the parameter errors 
have been rescaled by the reduced \/x^ value, are given 
in Table H 

With the present data it is difhcult to distinguish if 
this is the behavior of LhHz alone, or is due to blending 
with other components which are not strong and coher- 
ent enough to be observed separately on short time scales 
and are lost in the averaging of the power spectra. This, 
as well as other ambig uities (see discussion abo ut the 
identification of L^ow in Ivan Straaten et al.ll2003l ). arise 
because of the gaps between the u^, viow, Vh and VhHz 
versus relations (See Figure [S]). Although we are not 
able to explain those gaps, the interpretation that L/i^z 
could be affected by the presence of other components 
is made more likely by the fact that L{ in Z sources 



can be unambiguo usly identified down to frequencies of 
^ 150 Hz (see e.g. iJonker et al.ll2002bl) and that Lh and 
Liow have sometimes been observed at frequencies up 
to ~ 120 and ~ 60 Hz, respectively, i.e., in both cases 
reaching t he hectohertz Q PO range (van Str aaten et ahl 
[2005. Linares et al.l[2005[ ). If it would be possible to fol- 
low a source in its evolution from the extreme island 
states where L^ow is prominent, to the lower left banana 
where Lg is seen, then some of these ambiguities could 
be resolved. 





QhHz 




Slope 


(1.6 ±0.2) • 10-^ 


(-5.7 ± 1.1) ■ 10-^ 


Constant 


-0.52 ±0.15 


196 ± 10 


Amplitude 


0.35 ±0.04 


23.5 ±2.1 


Period 


532 ± 48 


498 ± 25 


X^dof 


15.1 / 15 


14.9 / 15 



TABLE 4 

The results of fitting the data of Qhhz and vhhz versus 
that of vu for the four atoll sources 4u 1728-34, 
4U 0614+09, 4U 1608-52 and 4U 1636-53. We used the 
combination of a straight line plus a sine. the parameter 
errors have been rescaled by the reduced y/y^ value (see 
text). The quoted errors use Ax'^ = 1.0. 




I 4U 1820-30 I • I i: 

4U 1728-34 I— . 
4U0614-^09 t 
4U 1636-52 
4U 1608-52 ^-m-^ 
0-1 - SAX J1 808.4-3658 " 

-.^ 1^ XTE J1 807-294^ I: 

200 400 600 800 1000 
(Hz) 

Fig. 10. — Llf's rms (above) and Q (below) versus Vu for 7 atoll 
sources. The symbols are labeled are the same in both plots. The 
scatter of the points are likely due to the averaging method (see 
text). 



4. DISCUSSION 

In this paper, we report a detailed study of the time 
variability of the atoll source 4U 1636-53 using RXTE 
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Fig. 11. — The fractional rms amplitude of all components (exc ept Li^p) plotted versus v ^. The symbols are labeled in the plot. The 
data for 4U 1728-34, 4U 1608-52 and 4U 0614+09 were taken from I van Straaten et all pOOSi '). Note that for L^^and of 4U 1608-52, 
the 3 triangles with vertical error bars which intersect the abscissa represent 95% confidence upper limits (see van Straaten et al. ( 200^ 
for a discussion). Also note that we exclude the 3 points in the rms£ versus Uu plot which were identified as Liow bv ivan Straaten et al] 
||200^) . 



that includes, for the first time, observations of this 
source in the (low-luminosity) island state. We divided 
the data into 15 intervals, A to N, based on the posi- 
tion of the source in the color diagram. Based on the 
fact that, (i) intervals Al, B, D and E show a narrow 
QPO with a characteristic frequency between i>h and i/h 
which was previously seen in other atoll sources when 
they were in their island st ate (e.g. Ivan Straaten et aTl 
l2002tlAltamirano et al.ll2005[ ): (ii) intervals Al, A2, B, D, 
E and F do not show either Lg or power-law VLFN at fre- 
quencies lower than 0.5 Hz, which would be expected to 



be present in the banana state (|Hasinger fc van der KlisI 
119891: Ivan der Klisll2004 |2006[) : (iii) the intensity of the 
source starts to increase from interval G (see Figure [2]) 
and (iv) intervals Al to F occupy the hardest loci in the 
color diagram (see Figures [5] and [2]) , we conclude that 
intervals Al to F show the source in the island state, 
representing the first RXTE observations of 4U 1636-53 
in this state. Interval G may represent the transition 
between the IS and the LLB since its power spectrum is 
very similar to that of the first five intervals, but with the 
difference that a weak (~ 1.2 % rms) VLFN is present 
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Fig. 12. — Quality factor Q of Lu , Li and Lf^fj^ versus i/u. 
Symbols are the same as in Figure 1111 Note that the results for 
the quality factor are probably affected by the averaging method 
(See Sections [T] [3] and the Appendix). 

at a frequency lower than 0.1 Hz (see Figured]). 

Along the color color diagram we find all seven power 
spectral components that wer e already seen in other 
sources in previous works (see Ivan der KlisI 12004 for a 
review): L„ is detected in all of our power spectra, Li 
is unambiguously detected starting from power spectrum 
H {vu ~ 800 Hz), LhHz is observed in 11 out of 15 power 
spectra at frequencies between 100 and 270 Hz, Lh and 
Llp are detected mainly in the island state, Lj, is always 
observed and finally Li,2 is detected when Lf, becomes 
peaked from interval G. 

Previous works have shown that the frequencies of the 
variability components observed in other atoll sources 
follow a unive rsal scheme of correlati ons when plot- 
ted versus i^u (|van Straaten et ahl I2003L and references 
therein). We have found that the noise and QPO fre- 
quencies of the time variability of 4U 1636-53 follow simi- 
lar correlations as well (see Section |3]) confirming the pre- 
dictive value of this universal scheme. However, we also 
found some differences between 4U 1636-53 and other 
atoll sources which we discuss below. As 4U 1636-53 
is one of the most luminous atoll sources showing the 
full complement of island and banana states (full atoll 
track), the object is of interest in order to investigate 



the luminosity dependence of spectral/timing state be- 
havior. This is of particular importance to the ongoing 
effort to understand the origin of the difference between 
the atoll sources and t he more luminous Z sources (see 
e.g. iHoman et aT1l2006l ). 
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Fig. 13. — 1^52 versus for the atoll source 4U 1636—53 (this 
paper) and the characteristic frequencies of the low-fr equency nois e 
(LFN) ve rsus Vu fo r the Z sources GX 17+2 (Homan et al.''2002l), 
Cyg X-2 IIKuznetso v 2002), GX 340+0 (Jonkcr ct al. 2000a)_and 
GX 5 -1 l lJonker et al.l I2002e]) — see text and Ivan Straaten et al.l 
pOO^ ) for details. (See also Figure IHll 



4.1. 



The broad components in 4U 1636-53 and Z- source 
LFN 

As can be clearly seen in Figure [U where we show 
V}j2 versus v^, the behavior of the L;,2 component dif- 
fers significantly between sources. For 4U 1636-53 and 
4U 0614+09, vi,2 increases with Vu, while this is not seen 
for 4U 1608-52 and 4U 1728-34. This frequency behavior 
is different from that observed for all other components 
(see Section 13]), which instead is very consistent between 
sources, even for the case of the hHz QPO, which has 
not been seen to correlate with other components (see 
Ivan der KlisI |2004| . for a review) . This unusual, some- 
what erratic behavior of L;,2 may be related to the fact 
that it is usually detected as a relatively weak wing to a 
much stronger L;,, so that small deviations in the time- 
averaged shape of L\j have a large effect on Lfc2- 

In order to investigate the relation of _Lf,2 to the well- 
known low frequency noise (LFN) which occurs in the 
same frequency range in Z sources, in Figure [T51 we plot 
the results for 4U 1636-53 toge ther with thos e for the 
LFN in th e Z sources (Hasinger fc van der K lis 198^ 
GX 11+2 (IHoman et alJ l200i. Gvg X-2 (iKuznetsovl 
l200l. GX 340-HO (|Jonker et all l2000al) and GX 5-1 
(jJonker et ani2002bD . Note that the broad-band noise 
in these Z sources was not fitted with a zero-centered 
Lorentzian but with a cutoff power law or a smooth 
broken power law. We used the results of the con- 
version from power laws t o zero -centered Lorentzians 
done bv Ivan Straaten et aP (l2003^. Previous works (e.g. 
iPsaltis et al.l 119991 : Ivan Straaten et al.l I2003D compared 
the time variability of Z sources with that of atoll sources 
and tried to associate variability components among 
these sources. Based on frequency-frequency plots, only 
the kHz QPOs and the horizontal branch oscillations 
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(HBO) found in the Z sources can be unambiguously 
iden tified with atoll sou r ce co mponents, the latter with 
Lh- Ivan Straaten et all (|2003f ) suggested that the LFN 
might be identified with Lf,2 and noted that (like in the 
case of Lb2) the characteristic frequency of the LFN, 
when plotted versus Vu-, does not follow exactly the same 
relations between Z sources. By comparing the differ- 
ent frequency patterns in Figure [131 we find that the 
behavior of the LFN component of GX 17+2, and that 
of Lf,2 of 4U 1636-53 are similar, which might indicate 
that the physical mechanism involved is the same. Per- 
haps this is related to the fact th at 4U 1636-53 is a rela- 
tively luminous atoll source (see lFord et ani2000[ ) while 
GX 17-1-2 may be a relatively low luminosity Z-source 
(|Homan et al.ll2006f ). Hence, 4U 1636-53 might be rela- 
tively close in to GX 17-f 2 and differ more in from 
the other two sources introduced above. Note that the 
time variability of GX 17-1-2 is different from that of the 
other Z sources plotted in Figure [131 For instance, the 
characteristic frequency of its LFN is rather low and it 
appears as a peak, it shows a flaring branch oscillation 
(FBO) and the harmonic of the HBO is relatively strong, 
whereas the other Z sources plotted show a flat LFN, no 
FBO and a weak harmonic to t he HBO (iJonker et al.l 
l2002bD . As previouslv noted bv iKuulkers et al.l (|1997D . 
these properties set GX 17+2 apart from the 'Cyg-like' Z 
sources GX 5-1, GX 340+0 and Cyg X-2 and, associate 
it with the 'Sco-hke' Z sources Sco X-1 and GX 349+2, 
not plotted in Figure [T3l because no systematic study of 
the of the LFN and QPO behavior of these sources in 
terms of Vmax is available as yet. 

We further investigated the frequency similarities be- 
tween 4U 1636-53 and GX 17+2 by plotting our results 
for the two sources. No clear component associations 
were found. GX 17+2 is the only Z source that had 
shown an anti-correlation between the frequency of one 
of its compon ents (HBO) and th e kHz QPOs at high 
Vu > 1050 Hz (jHoman et al.|[200^ . A similar effect was 
seen in the atoll source 4 U 0614+09 between i/f, and Vu 
(Ivan Straaten et"ani2002f ). As can be seen in Figure \5\ 
a similar decrease of with z/„ at high frequency may 
occur in 4U 1636-53. However, the error bars on Vb are 
rather large in the relevant range, and the data are still 
consistent with vi, being constant at ^ 1100 Hz, and 
marginally, even with a further increase in frequency. It 
is interesting to note, that while 4U 0614+09 has a much 
lower Lx than 4U 1636-53, both sources might show this 
same turnover in Vf, versus Vu- Of course, these results 
need confirmation. 

4.2. The low frequency QPO 

With respect to the low-frequency QPO L^p, 
Ivan Straaten et all ()2003L [2005 ) observed that in their 
data there were two groups of sources, one where the 
Llf feature was visible, and a second one, were a 
QPO was detected which they suggested to be the sub- 
harmonic of Llf and therefore, called Following 
Ivan Straaten et ahl (|2003f l. the upper continuous line in 
Figure [9l indicates a power law fitted to the i^lf versus 
Vh relation of the low luminosity bursters IE 1724-3045 
and GS 1826-24, and the BHC GX 339-4. If we repro- 
duce the fit where we take into account the errors in both 
axes, we find a best fit power-law index a ~ 0.97 ± 0.01 
and x'^/dof = 80/19 - 4.2. If we fix a = 1, the fit gives 



a X^/dof = 83/20 ~ 4.1. According to the F-test for 
additional terms, there is a < Ict improvement of the fit 
when a is set free, so we conclude that vlf is consistent 
with being linearly related to Vh- The lower dashed line 
is a power law with the same index a = 0.97, but with a 
normalization half of that of the dashed line. 

As can be seen in Figure [SI in 4U 1636-53 the L^p 
component does not follow either of the two power-law 
fits. If we fit the points for 4U 1636-53, we find that 
the power-law index is a2 = 1-40 ± 0.09 {x^ /dof — 
0.14/2), significantly different from that of the other 
sources. Given the above, it is probably incorrect to 
think that the difference in the vlf vs. vt relation be- 
tween GX 339-4, GS 1826-24 and 1724-3045 on one 
hand and 4U 1608-52, Cyg X-1 and XTE J0929-314 on 
the other is associated w ith harmonic mode switching 
(jvan Straaten et aLll2003f). This co nclusion is supported 
bv the work of iLinares et al] (|2005f ) who also found a dif- 
ferent correlation [a = 0.58 ± 0.06, see also Figure [U in 
XTE 1807-294 over a much wider range of frequencies 
than we obtained for 4U 1636-53, by the high x^ /dof 
for the Vlf — fit on the data of the low luminos- 
ity bursters IE 1724-3045, GS 1826-24, and the BHC 
GX 339-4 (see previous paragraph), by the fact that if 
we use the centroid frequencies instead of Vmax, the re- 
lations worsen (see j/an Straaten et a l. 2003i'l, and by the 
fact that the points for 4U 1728-38 (jvan Straaten et aTl 
120021 ) fall in between the two power laws, (solid and 
dashed line in Figure [9]) . Nevertheless it is interesting 
that the data for 4U 0614+09, 4U 1728-34, 4U 1636-53, 
4U 1820-30, 4U 1608-52 XTE 1807-294, SAX J1808.4- 
3658, XTE J1814-338 and XTE J0929-314 all fall on, or 
in between, the two previously defined power laws, i.e., 
do not deviate from a single relation by more than a fac- 
tor of 2. We note that all the v^p values discussed here 
could in principle have been affected by smearing in the 
averaging process discussed in Section [2l However, for 
smearing to shift a frequency-frequency point away from 
its proper value, large systematic differences are required 
between the two components in the dependence of am- 
plitude or Q on frequency, and in the case of L^p and 
Lh there is no evidence for this. 

From Figure [9l it is apparent that the frequency 
range in which the L^p component has been iden- 
tified is rather large (up to 2.5 decades). Clearly, 
which frequency ranges L^p covers is not related to 
source spin frequency, angular momentum or luminos- 
ity of the object. The sources 4U 1608-52, 4U 1820-30, 
4U 1636-53 and 4U 1728-34 all show Llf when they 
are in their island state, but with v^p < 2.6 Hz for 
4U 1608-52, and > 30 Hz for the other three sources. 
The accreting millisecond pulsar XTE J1807-294 shows 
VLF > f2 Hz while the AMP XTE J0929-314 shows 
i^LF ^ 1 Hz, wh ile both have very sirn ilar spin fre- 
quencies fl91 H z. [iviarkwardt et al.l l2003l and 185 Hz, 
iRemillardet al.l |2002| . respectivelvl . 4U 1820-30 and 
4U 1636-53 are at least one order of magnitude more 
luminous than XT E 1807-294 and SAX J1808.4-36 58 at 
their brightest fsee lFord et al.ll2000l:lwiinanddl2005l ). but 
all four sources show vlp > 5 Hz. The only systematic 
feature in the LF QPO frequencies is that while frequen- 
cies up to 50 Hz are seen in neutron stars, black holes 
have not been reported to exceed 3.2 Hz, nor did atoll 
sources in the EIS exceed 2.6 Hz. So, BHCs and NS 
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in the extreme island state are similar in this respect; 
(this may be related to an overall similarity in power 
spectral shape for such sources in these stat es that was 
noted before; see, e.g., [P saltis et al. 1999; Nowaldl200(il : 
iBelloni et aIll2002Hvan Straaten et al.ll200l . 

4.3. The X-ray luminosity dependence of rms 

It has been suggested that an anti-correlation may 
exist between the average X-ray luminosity of differ- 
ent sources and the rms amplitude of their power spec- 
tral co mponents (see discussion in iJonker et al.l 120011 : 
Ivan St raaten et al. 2002, 2003, and references therein). 
From Figure [TTj wc find differences in kHz QPO rms am- 
plitudes of no more than a factor 2 between sources which 
differ in average luminosity by a factor up to 10, except 
for one point of 4U 0614-1-09 at Vu ^ 1140 Hz, where the 
rms of the upper kHz QPOs is a factor ~ 7 h igher than 
that of the other atoll sources. iMendez et al. (2001), 
IJonker et all ()2001f ) and Ivan Straaten et al.l (|2002,) have 
already noted that the data are inconsistent with a model 
in which the absolute amplitudes of the kHz QPOs are 
the same among sources, and the decrease in rms with lu- 
minosity between sources is only caused by an additional 
source of X-rays unrelated to the kHz QPOs. 

From Figure [TT] it can also be seen that the largest 
rms amplitude differences are found in the hHz QPOs 
(excluding 4U 1608-52, which is a transient source cov- 
ering a large L,j: range). For this component we find (1) 
4U 0614-f09, which has the strongest rmshHz (> 15% 
when Vu < 800 Hz); (2) 4U 1636-53, which has the 
weakest rmshUz (< 10% when Uu < 800 Hz) and (3) 
4U 1728-34 which has rmstHz generally between those 
of (1) and (2) [between 10 and 15% when < 800 Hz]. 
(At Vu > 800 Hz, the groups can still be differentiated as 
the rms amplitud e decreases with Vy)). From figure 1 in 
iFord et all pOOOl ). it can be seen that 4U 0614-^09 is the 
faintest X-ray source of our sample, while 4U 1636-53 
is the brightest. 4U 1728-34 show luminosities between 
the first two. This suggests an X-ray luminosity-rms 
anti-correlation for LhHz that is not as clear in the other 
components (see also Figure [TU|) . 

The fact that the rms of LhHz starts to decrease at the 
same as that of Lu and Lb^ while VhHz does not corre- 
late with i^u as all other frequencies do, suggests that the 
frequency setting mechanism is different for LhHz com- 
pared with the other components, while the amplitude 
setting mechanism is common. As pointed out in Sec- 
tion [31 the drop in rms in Lu, LhHz and L^ starts at z/„ 
between 700 and 800 Hz. For the case of 4U 1636-53 
shown here, this corresponds to interval G. The power 
spectrum of this interval may represent the transition 
between the island and the banana state, when the geo- 
metr ic configur ation of the syst em is thou ght to change 
(e.g. IJonker et al. 2000b; Gierli hski fc Do ne 2002). For 
example, the appearance of a puffed-up disk could smear 
out the variability coming from the inner region where 
the oscillations are produced. 

4.4. The nature of the hectohertz QPOs 

While our results indicate that the characteristic fre- 
quency of the hHz QPO may oscillate as a function 
of VhHz remains constrained to a limited range 
of frequencies (100-250 Hz) for 4U 1636-53 and for 



the other sources used in Figure [5l A similar re- 
sult has been reported for VhHz i n several other ato ll 
sources such as in MXB 1730-335 dMighari et al.l l2005D. 
4U 1820-30 (jAltamirano et al.l [200'5f ) and in the atoll 
sourc e and millisecond accr eting pulsar SAX J1 808.4- 
3658 (jWiinands fc van der Klis 199 8: van Straaten et al.l 
I2005D . Interestingly, the presence of LhHz bas not been 
confirmed for Z-sources ()van der Klisll2006l ). possibly due 
to the intrinsic differences between atoll and Z-sources 
s uch as luminosity. 

Ivan Straaten et ah! (|2002|) have sugges ted a link be - 
tween the < 100 Hz QPOs reported by iNowakI ()2000[ ) 
in the black holes C yg X- 1 and GX 339-34 and LhHz- 
Ivan Straaten et all (|2002t ) also suggested that LhHz 
could be related to the ~ 67 H z QPO in the black hole 
GRS 1915-1-105 (IMorgan et al.lir997l and the ^ 300 Hz 
QPO in the BHC GRO J1655-40 (tlemillard et al.l | 1999|) 



20011) 



which also have stable frequencies. iFragile et al.T 
made a tentative identifi cation of the ^ 9 Hz QPO in 
the BHC GRO J1655-40 (Remill ard eFaLlflggg! ) with the 
orbital freque ncy at the Bardeen-Petters on (B-P) tran- 
sition radius (|Bardeen fc Pettersonlll975l ) and suggested 
the same identification for LhHz in neutron star systems. 
In this scenario, the orbital frequency at the radius where 
a warped disk is forced to the equatorial plane by the 
Bardeen-P etterson effect can produce a quasi-periodic 
signal (see IFragile et al.]|2001l for an schematic illustra- 
tion of the scenario). 

Attempts have been made to theoretically esti- 
mate the B-P transition radius from accretion disks 
models in terms of the angular momentum and the 
mass of the compact object (e^ . IBardeen fc PettersonI 
19751 llvanov fc Illarionovl 119971: [Hatchett et all 119811: 



Nelson fc Papaloizoul |2000|) . IFragile et all ()200lD pro- 



pose a parameterization involving a scaling parameter 
A, which according to them lies in the range 10 < A < 
300. These authors write the B-P radius as Rbp = 
A ■ a^J'^ ■ Rgr, where a* = Jc/GM^ is the dimension- 
less specific angular momentum (J and M are the an- 
gular momentum and the mass the compact object, re- 
spectively) and Rgb. is GM/(?. The Keplerian orbital 
frequency associated with the B-P transition radius can 
be written as Vkep,BP = ■ (27rG)-i • {Ma^^A^^^y^ . If 
we assume that the atoll sources plotted in Figure [5] all 
have m asses between 1.4 and 2Mp-,, that 0.3 < a* < 0.7 
(see e.g. lSaTgado et al.|[l99llCook et al.lll994 and refer- 
ences within) and that the central frequency of the hHz 
QPOs is between ~ 100 and ~ 250 Hz, we can con- 
strain the scaling factor A for these source to be between 
20 and ^ 84. If A only depends on the accretion 
disk (i.e. does not depend on the central object), this 
can be used to constrain the frequency range in which 
we expect to observe Vkep,BP in black holes. For ex- 
ample, the black hole BHC GRO J1655-40, whose mass 
is estimated from optical and in frared investigations as 
M = 6.3±0.5Mq ([Greene et al.ll2001. ') and whose specific 
angular mom entum can be estimated to be between 
0.5 and 0.95 (|Cui et al.lll99l: [Fragile et al.ll200l . would 
have i'kep,BP between ~ 6.5 and ~ 127 Hz, which would 
exclude the 300 Hz QPO observed in GRO J1655-40 but 
wo uld be consistent wi th the the 9 Hz QPO as proposed 
bv IFragile et"al1 (|200]1 ). If one assumes the 450 Hz QPO 
in GRO J1655-40 is associated with orbital motion at 
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the last stable orbi t, then a^, could be as low as ^ 0.15 
()Strohmaveill2001bh . In this case, Vkep,BP can be as high 
as ^ 425 Hz for a black hole mass of 5.7M0. 

5. SUMMARY 

• Our observations of 4U 1636-53, including the first 
RXTE island state data of the source, show timing 
behavior remarkably similar to that seen in other 
atoll NS-LMXBs. We observe all components pre- 
viously identified in those sources, and find their 
frequencies to follow similar relations to those pre- 
viously observed. This is interesting as the sources 
compared in this work were observed at intrinsic lu- 
minosities different by more than an order of mag- 
nitude. 

• The previously proposed interpretation of the QPO 
frequencies v^p and i^lf2 in different sources in 
terms of harmonic mode switching is not supported 
by our data on 4U 1636-53, nor by data previously 
reported for other sources. However, these frequen- 
cies still do not deviate from a single relation by 
more than a factor of 2 for all sources. 

• The low frequency QPO Xlf is seen in black holes 
and in accreting millisecond pulsars as well as in 
non-pulsing neutron stars at frequencies between 
^0.1 and ~ 50 Hz. The frequency range that 
Llp covers in a given source is not related to spin 
frequency, angular momentum or luminosity of the 
object. 

• The rms and frequency behavior of the hectohertz 
QPO suggests that the mechanism that sets its fre- 
quency differs from that for the other components, 
while the amplitude setting mechanism is common. 
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Appendix 

In this Appendix we further discuss other possible ap- 
proaches to analyze the characteristics of complex power 
spectra such as generally found in neutron star low-mass 
X-ray binaries. 

In the ideal case, we would have data with enough 
statistics to be able to follow the evolution of the parame- 
ters of all observable components in the power spectra on 
sufficiently short time scales to be sensitive to the small- 
est meaningful variations. Unfortunately, this is not the 
case for the present data and meaningful variations are 
averaged out in our data. These variations can some- 
times be recovered by the use of alternative methods. 
For exa mple, with the "shift a nd add" method intro- 
duced by lMendez et al.l (|1998bl ). it has been possible to 
better constrain some of the characteristics of the kHz 
QPOs in several sources than without this method (e.g. 



iMendez et allllQQSbl : iBarret et al.ll2005H) . A disadvan- 
tage of the method is that it distorts the power spectrum 
at the lowest and highest frequencies covered. 

We investigated if this method could also be used for 
our purpose. However, in our experiments with this we 
encountered several complications. From the observa- 
tional point of view, in order to use this method we re- 
quire a sharp power spectral feature that can be accu- 
rately traced in time. There are two possibilities for such 
features: the lower kHz QPO and the low-frequency 
Lorentzian LiF. While the lower kHz QPO is usually su- 
perimposed on well- modeled Poisson noise, tracing L^p 
is complicated by the fact that it i s superim posed on 
strong variable broad band noise (see Ivan der Klis 200(|, 
and references within). More importantly, while Li can 
be traced on sufhciently short time scales (< 64sec) for 
the intrinsic changes in the characteristics of the power 
spectrum to be minimal, typically an entire observation 
is required for detecting L^p at sufficient signal to noise. 
In practice this means that we can only use the shift and 
add method with L^, which constrains us to only that 
relatively limited part of the data where Li is actually 
detected (see Figure [5]). We note that Lpp and are 
not simultaneously detected in our data set. 

We analyzed all the datasets described in Section [2] 
and found that ^ 15% (^ 0.17 Msec) of our data have 
traceable lower kHz QPOs. Most of that time the lower 
kHz QPOs are detected at frequencies between 700 and 
850 Hz (the full range was 600-900 Hz). 

In order to use the shift and add method we must adopt 
a relation between the frequency of the component we 
wish to shift on (here the lower kHz QPO) and the fre- 
quency of the component we wish to detect (here the low 
frequ ency QPOs/noise). F or example, in their original 
work lMendez et ahl (|1998bD supposed that the difference 
between the lower and the upper kHz QPO frequency 
remained constant when both peaks move. The study 
of the characteristics of the low frequency QPOs/noise 
using the shift and add method is complicated by the 
fact that we have imprecise information about their re- 
lation with the lower kHz QPOs: a constant frequency 
difference certainly does not apply even to narrow ranges 
in shift frequency. The aim of this paper as well as the 
aim of the papers cited below is to present observational 
result s that help constra i n thos e relations. 

As Ivan Straaten et al.l (j 20021 [2003) showed, the fre- 
quencies of all components except those of the hHz 
QPOs are correlated in a similar way betw e en so urces 
(see Figure |5] ). How ever, [ van Straaten et all (j2005f ) and 
iLinares et al.l (|2005l ) also showed that those correlations 
are sh ifted in pulsating sources and lAltamirano et ahl 
(|2005f ) found that even non-pulsating systems might 
show frequency sh i fts. Additi o nally, t he results of 
lyaii_Straaten et all (12001 l2003l. [2005( 1 . ILinares et all 
(|2005D and lAltamirano et al] (|2005l ) show that although 
the frequency relations between the different components 
are well fitted with a power law, the index and normal- 
ization of the power law are different for each relation 
and may also depend on frequency range. 

In order to quantify the problem described above, we 
studied observation 60032-01-05-00 using power spectra 
of 64-sec data segments at 2 Hz frequency resolution. 
This is a very good observation for our purpose since: (i) 
it has ^ 27 ksec of uninterrupted data; (ii) the lower kHz 
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QPO is strong enough to be significantly detected within 
64 seconds for the entire 27 ksec; (iii) the lower kHz QPO 
frequency drifts between ^ 700 and ~ 860 Hz and (iv) 
the power spectrum can be fitted with 4 Lorentzians: 2 
for the kHz QPOs, one for if, (at 40 ± 4 Hz) and one for 
Lb2 (at 20± 12 Hz; 2.9a). 

We first analyzed the power spectrum obtained by 
aligning the Li components. We found that Lb and Lh2 
had blended into a broad component at ~ 100 Hz. This 
result was expected, as a drift of 160 Hz in vi does not 
imply a drift of the same magnitude in the frequencies of 
the low v components. We then tried to align the power 
spectra by predicting the position of the low-frequency 
components from vi using a different power law relation 
for each component as reported by Ivan Straaten et al.l 
([2005^ for Lb and Vb2 = e-24±8 . j,+3.8±i.3 ^j. (^j^g j.g_ 

lation for Lb2 is based on our data for 4U 1636-53; given 
the large errors in our data, the j doj for the power law 
fit was 0.26) . The results depended on which power law 
we used: no significant changes in the resulting power 
spectrum were found when trying to align Lbi , while the 
power of both components was smeared out producing a 
blend when we tried to align Lb- This was clearly the 
effect of the difference in power law indices and normal- 
izations between components. As the frequency relations 
we used are between and the frequencies of the other 
components, we had to assume a relation between vg and 
in order to predict the frequency variations. We var- 
iously assumed l\v = — ^ 280, 300 & 350 Hz and 
obtained similar results in each case. We also predicted 
by fitting a line to the Az^ vs. data reported by 
iJonker eT al.r(|2002af) in the range 720 < vt < 900 Hz. 
Again, the results of the power-spectral fits were the 
same within errors as those of the previous experiments. 

We repeated the last exercise (using the power law re- 
lations) also for all 0.17 Msec of vi useful data and for 
Lf,2, Lb a nd Lb (we use the power law relation as re- 
ported bv Ivan Straaten et al.l l2005l for L^). We again 
found that our results were dependent on the power law 
used and not significantly better defined than the aver- 
age power spectra obtained when all 0.17 Msec of data 
were averaged together without shift. 

In another experiment we calculated 4 average power 
spectra including all 0.17 Msec of useful data by select- 
ing only those 64-second segments which had between 
650-700, 700-750, 750-800, and 800-850 Hz respectively, 
and averaging these selected power spectra without shift- 
ing. In all cases we detect both kHz QPOs, a power 
law VLFN and Lb. As expected from the results shown 
in Figure |5l Vb is correlated with the frequency of both 
kHz QPOs. The measured frequencies are all consistent 
within errors with those reported on Figure (5) and Ta- 
ble [5] The lack of statistics in each average power spec- 
trum did not allow us to well constrain the power spectral 
parameters of other components. 

Finally, we fitted a line to the relation between vj^p 
and Vh defined by the 4-points visible in Figure |9l We 
used the vj^p we find in all four power spectra (Al, B, D 
and E) to predict and shift and add these four power 
spectra together. Again, we find that the blend of com- 
ponents (this time between Lb and L^) and the distortion 
of the power spectra at low frequencies prevented us to 
better estimate the L^ parameters. 

So, neither the shift and add method nor selecting data 



on in 64-sec segments (i.e. much shorter than an obser- 
vation) in our data provides an advantage in measuring 
the broad low frequency components better. Therefore, 
in this paper we decided to the straightforward method 
described in Section |2l 
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TABLE 5 Characteristic frequencies Umax, Q values 

(fmax = fcentral/FWHM), FRACTIONAL RMS (iN THE FULL PCA ENERGY 
band) and component IDENTIFICATION (ID) OF THE LORENTZIANS FITTED 

FOR 4U 1636-53. The quoted errors use Ax^ = 1.0. Where only one 

ERROR IS QUOTED, IT IS THE STRAIGHT AVERAGE BETWEEN THE POSITIVE AND 
THE NEGATIVE ERROR. NOTE THAT THE RESULTS FOR THE QUALITY FACTOR OF 
BOTH AND Li^p COMPONENTS ARE AFFECTED BY OUR AVERAGING METHOD 

(See Sections [T] [3] AND the Appendix) 



l^max (Hz) 


Interval Al 

Q RMS (%) 


ID 


434.8 ±26.1 

175.9 ±8.9 
15.0 ±0.6 

5.1 ± 0.1 

3.2 ±0.4 


1.6 ± 0.5 

2.7 ± 1.3 
0.7 ±0.1 

7 q+9.0 
— 2 2 

0.33 ± 6.07 


12.9 ± 1.4 
8.06 ± 1.23 
13.3 ±0.7 
3.2 ± 0.6 
11.1 ±0.7 


LhHz 

Lh 
L, 


I'max (Hz) 


Interval A2 

Q RMS (%) 


ID 


436.4 ± 28.2 
12.3 ±0.8 
2.34 ±0.52 


0.92 ± 0.24 
0.27 ±0.11 
0.14 ±0.09 


16.5 ± 1.05 
16.1 ± 1.1 
9.9 ± 1.4 


Ll 
Lh 


J^max (Hz) 


Interval B 
Q RMS (%) 


ID 


464.9 ±9.4 

18.2 ±0.3 
6.83 ±0.11 
4.03 ±0.31 


1.7 ±0.2 
0.53 ± 0.42 
0.75 ± 0.06 
3.07 ±0.66 
0.19 ± 0.03 


13.8 ± 0.9 
7.4 ± 1.8 
12.4 ±0.4 
4.1 ± 0.5 

11.9 ±0.4 


LhHz 

Lh 
Llf 
Lh 


i^Tnax V / 


Interval C 
Q RMS (%) 


ID 


529.3 ± 15.4 
23.1 ± 1.6 
6.38 ± 1.01 


1.1 ±0.1 
0.47 ±0.15 
0.09 ± 0.05 


17.1 ±0.6 
11.9 ± 1.5 
13.1 ± 1.2 


Lu 

l\ 


i^max 1 


Interval D 
Q RMS (%) 


ID 


524.7 ±8.0 
201 7+''8-0 

23.9 ± 0.5 
9.8 ±0.3 
5.6 ±0.7 


2.2 ± 0.4 
0.58 ± 0.35 
0.99 ± 0.13 
2.17 ±0.63 
0.15 ± 0.04 


13.9 ± 1.0 
8.6 ± 1.9 
10.5 ±0.6 
5.2 ± 0.9 
11.4 ±0.7 


-^it 

LhHz 

Lh 
Lh 


l^max (Hz) 


Interval E 
Q RMS (%) 


ID 


593.1 ±4.3 
270.3 ± 31.7 

31.2 ± 1.3 
15.04 ±0.41 

8.5 ±0.9 


3.1 ±0.2 
0.63 ± 0.20 
0.89 ± 0.15 

1.6 ±0.5 
0.22 ± 0.03 


12.7 ±0.4 
8.7 ±0.8 
8.7 ±0.7 
5.7 ± 1.2 
11.2 ±0.6 


LhHz 

Lh 

•^LF 
Lh 


t^max (Hz) 


Interval F 
Q RMS (%) 


ID 


637.2 ± 7.9 
208 7 ± 30 1 
18.7 ±0.9 


3.06 ±0.32 

1.2 ± 0.6 
0.02 ± 0.04 


15.9 ±0.6 
7.2 ± 1.2 
16.9 ±0.2 


^hH z 

Lh 


l^max (Hz) 


Interval G 
Q RMS (%) 


ID 


78U.01 ± 1.78 
152.7 ±10.4 
23.4 ±0.5 
3.8 ±0.9 


4.5 ± 0.1 
0.32 ± 0.11 
0.34 ±0.03 
0.20 ±0.11 


13.6 ±0.1 

8.7 ±0.49 
11.6 ±0.3 

2.8 ± 0.5 


LhHz 
Lb 


I'max (Hz) 


Interval H 
Q RMS (%) 


ID 


862.8 ± 1.& 
568.3 ±4.6 
120.6 ±3.6 
26.6 ±0.5 
4.3 ±0.7 


6.2 ±0.2 
7.06 ± 1.22 
1.03 ± 0.14 
0.57 ±0.04 
0.35 ± 0.11 


11.4 ±0.1 

4.6 ± 0.2 
6.4 ±0.3 
9.2 ± 0.2 

2.7 ±0.3 


Lu 
Li 

l^hHz 

Lb 


l^max (Hz) 


Interval I 
Q RMS (%) 


ID 



585.5 ± 5.3 
228.81 ± 9.05 
112.68 ± 4.34 
28.6 ±0.6 
6.06 ± 1.41 



7.6 ± 1.1 

5 6+5-6 

1.51 ± 0.33 
0.80 ± 0.09 
0.32 ± 0.14 



5.3 ± 0.3 
2.6 ±0.5 
5.29 ± 0.41 
7.72 ± 0.28 
2.93 ± 0.41 



Le 

J harmonic 
^hHz 
LhHz 





Interval K 






r\,ivi o ) 


ID 


998.18 ± 7.94 
728.04 ± 2.63 
148.2 ± 8.1 
38.8 ± 0.9 


5.76 ±0.74 6.14 ±0.23 
5.3 ±0.3 8.3 ±0.2 
1.6 ±0.5 3.4 ±0.3 
1.2 ±0.2 4.69 ±0.37 


Le 




0.07 ±0.35 2.48 ±0.79 


Lb2 




Interval L 




y-max (Hz) 


Q RMS (%) 


ID 



1138.37 ±9.02 
837.5 ± 1.0 
156.1 ± 33.6 
46.4 ±0.9 
14.9 ± 5.8 



-L-±.VJO_3 gg 

11.4 ±0.4 

1.3 ±0.5 

2.4 ±0.4 
0.00 ±0.00 



2.7 ±0.3 
7.9 ± 0.09 
2.03 ± 0.38 

2.8 ± 0.2 

1.9 ± 0.3 



I'max (Hz) 

1220.5 ± 11.2 
44.2 ± 7.2 



Interval M 
Q RMS 

16.6 ±6.5 T7T 

0.56 ±0.28 3.1 



0.3 



Interval N 
(Hz) Q RMS (%) 



1259.1 ±9.9 
36.7 ±4.9 



^+a4.V 
-5.4 

0.05 ±0.13 



14.4"* 



1.8 ± 0.2 
2.37 ±0.11 





Interval J 




I'max (Hz) 


Q 


RMS (%) 


ID 


971.1 ±4.7 

666.2 ± 2.1 
127.9 ±9.6 
34.6 ±0.8 
8.07 ± 1.56 


9.9 ±1.1 
8.7 ±0.5 
1.1 ± 0.3 
1.1 ± 0.1 
0.34 ±0.14 


6.7 ±0.2 
7.7 ±0.2 
4.3 ± 0.4 
5.5 ± 0.2 
2.99 ± 0.31 


L, 

J^hHz 
Lb 



l^hHz 

Lb 

J±b2. 



ID 



ID 



( ;nntiniind on next ( :nliimn~l 
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